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ABSTRACT 


A study and pursuant development is described of a digital computer 
program for the computation of transfer function sensitivities in a 
symbolic form suitable for storage and subsequent repetitive numerical 
evaluation over a range of frequencies or parameter values. The algo- 
rithm is implemented in conjunction with the Network Analysis for Systems 
Applications Program (NASAP) version developed by R. S. Schwartz at 
Northeastern University. Several example problems, as well as sugges— 
tions for possible improvement of the algorithm and its extrapolation 


into the areas of optimization and automatic design, are included. 
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I. INTRODUCTION 


Pye STATEMENT OF THE PROBLEM 

The purpose of this research was twofold: (1) to investigate the 
feasibility of providing for the automatic digital machine computation 
of transfer function sensitivities in a symbolic form suitable for 
storage and subsequent numerical evaluation over a range of frequencies 
and parameter values or for use by a minimization/optimization procedure 
as a readily available source of gradient data, and (2) to implement 
such a program, if feasible, in conjunction with the Network Analysis for 
Systems Applications Program (NASAP). 

There were formidable obstacles encountered in the course of the 
research which may best be illustrated by a preliminary examination of 
the computer core storage required for implementation of the proposed 
procedure if it were meshed with an existent circuit analysis program 
such as NASAP. The NASAP program is ea pani of manipulating up to 10,000 
independent loop products (generated through use of topological flowgraph 
techniques), once as a AumeBecer polynomial in s and once as a denominator 
polynomial; that is, by iterative use of a single vector variable, LOOP, 
dimensioned to 10,000. Each loop product may represent from one to 
twenty-three constituent elements, and the analytical expression for the 
partial derivative of a quotient involves four separate components - 
numerator polynomial, denominator polynomial, partial derivative of the 
numerator and partial derivative of the denominator. The implication, then, 


is a storage allocation requirement of four vector variables each of 





dimension 230,000 words (920K bytes). This, of course, would provide 

for the sensitivity with respect to a single parameter; if the gradient 
is desired the requirement for the two partial derivative variables 
would be multiplied by a factor of 23 again, for a total of 11,500,000 
words. Consider now an even more staggering statistic. Since the NASAP 
program, like most programs written in the FORTRAN language, assigns 

only one alphanumeric character to a word (four bytes in the IBM 360/67), 
there would be an additional factor of four; three for the element 
identifying characters and one for a code to indicate the algebraic 
formulation required. This initial examination, then, would indicate the 
need for an impossibly high storage allocation of 46,000,000 words 
(184,000K bytes!). Some degree of success was achieved in the reduction 
of this figure to a manageable size, but at the expense of limitations 
on the original objective. The extent of the progress and associated 


restrictions will become clear in later sections. 


B. DISCUSSION OF THE NASAP PROGRAM 

Prior to discussing the sensitivity algorithm in any detail, it 
would perhaps benefit the reader to summarize the characteristics and 
capabilities of NASAP since it is upon this foundation that the proposed 
technique is firmly based. 

The revised version of NASAP is a digital computer program for the 
flowgraph analysis of lumped, linear active electrical networks designed 
to provide any desired circuit transfer function in either symbolic or 
numerical form. The program utilizes the "mixed" method to generate 
circuit equations, i. e., operations are performed on a tree-branch, 


link-branch topology so that Kirchoff's voltage law may be utilized to 





oa “ae «oe ao em om 


— 





express each link voltage in terms of tree voltages and Kirchoff's 

current law may be used to express each tree current in terms of link 
currents. These relationships, in association with volt-ampere and source 
dependency relationships, comprise a formulation defining a canonical form 
called a primitive flowgraph [3]. 

The user of NASAP prepares a circuit for analysis by: (1) converting 
to a linear model if necessary, (2) labeling each circuit element by type 
(R,L,C,V,I) and a two digit number which uniquely identifies it, and (3) 
numbering. each node as desired except for the ground node which must be 
designated node 1. The user then inputs this connectivity data in the 
following format: 


ELEMENT FROM TO ELEMENT CONTROLLING PARAMETER 
NAME NODE NODE VALUE (For controlled source) 


Following this are coded sets of data describing the transfer functions 
desired, in the form of the ratio of a response variable (voltage or 
current associated with some circuit element) to an excitation variable 
(voltage or current), and the forms of program output desired. The NASAP 
tree-~link hierarchy [20] and the order in which the elements were input 
to the program determines which elements are placed in the tree, thus 
providing the user with some flexibility regarding its choice.-. Given the 
topology, the machine then forms the fundamental outset matrix and from 
this constructs the primitive flowgraph. 

A transfer function between any two nodes of the primitive flowgraph 
may be obtained by adding a "closure" [20] and solving the resulting 
closed" graph by application of Mason's gain formula in the following 
manner. If the transfer function between node j and node k of a flow- 


graph is equal to ie and if the graph is then closed by the addition of 
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a branch from k to j whose transmission gain is equal to the reciprocal 
of Tay: then the loop gain with the added branch is equal to unity and 
the determinant of the flowgraph vanishes. This reduces the solution of 
7 to the problem of iin loop-set transmittance products, a procedure 
highly suitable to machine solution. Readers desiring more detailed 
information concerning this procedure are directed to Ref. 20. 

The procedures utilized by NASAP for labeling, sorting and storing 

each loop-set consist primarily of: 

(1) As each loop-set is determined by the search routine, a signed 
integer label is generated from which the loop-set may be 
uniquely regenerated. 

(2) Loop-sets with the same label but opposite sign are eliminated 
to reduce storage requirements. 

(3) Loop-sets are sotted and stored in order of decreasing magnitude. 

(4) The completed loop-set array represents the symbolic solution of 
either the numerator or denominator oo lnchetanen coded form. 
From successive solutions the transfer function may be evaluated 
either numerically or symbolically by decoding and summing each 


loop. 


Cx A COMPARISON OF TWO FAVORABLE TOPOLOGICAL METHODS 

A topological method other than the flowgraph procedure just described 
is the so-called "k-tree" approach [6] which is based on Maxwell's topolo- 
gical formulas for node voltage equations. That is, transfer functions are 
derived from the evaluation of determinants formed from circuit nodal 
equations;however, they may equivalently be formed in a topological fashion 
by summing all tree admittance products and appropriate 2-tree, 3-tree, ... 


k-tree admittance products. 





In both methods, k-tree and loop-set, the solution of passive 
circuits is obtained as a sum of linear products of the network parameters, 
thus eliminating the requirement for operations with polynomials, and 
facilitating machine evaluation of symbolic solutions. In addition, the 
two methods compare favorably with respect to the number of loop-trans- 
mittance products and tree-admittance products for a given circuit, since 
the maximum number of trees possible for a given topology is only one 
greater than the maximum number of loop-sets contained in the corresponding 
primitive flowgraph [20]. Each tree-admittance product is unique; however, 
in the flowgraph method several loop-sets are frequently found to have 
identical transmittance product magnitudes, although only one will remain 
after cancellation. Thus each remaining product is unique. A primary 
advantage of both procedures is this feature which obviates the evaluation 
of zero terms, as they would be if the system determinant were evaluated 
by matrix algebra. 

Although the two methods appear equal from the aieve discussion, a 
very important difference arises in the solution of active networks since 
the k-tree approach requires the graph to be modified and divided into 
separate current and voltage graphs from each of which the complete k-trees 
must determined, whereas in the loop-set analysis the procedure is invariant 
for passive or active networks. 

Another significant disadvantage of the k-tree analysis is its 
dependence upon the particular transfer function. to be evaluated, i-.e., 
impedance, admittance, gain, etc., each requiring different types of k-tree 
searches. Thus, the number of k-tree searches will vary with the number of 


transfer functions as well as the type desired. Therefore, the flowgraph 





method appears to have a distinct advantage in this respect, since the 
procedure remains identical regardless of the number or types of transfer 
functions. For these reasons the flowgraph approach was chosen for imple- 


mentation in NASAP. 


DD. SUMMARY 

It was the objective of this study to determine the feasibility of 
utilizing the loop-set symbolic form of the transfer function originated 
in the manner described in (B) above as a basis for computing and storing 
sensitivity information in similar symbolic form, thereby vastly simplify- 
ing the repetitive computation of accurate network sensitivities and 
rendering such solutions free of the difficulties involved with differen- 


tial, or numerical techniques. 





II. DEVELOPMENTAL PROCEDURE 


Utilizing the revised version of NASAP developed for the PDP-10 and 
CDC-3300 computers by R. S. Schwartz [20], adaptation to the IBM 360/67 
computer at the Naval Postgraduate School, Monterey was accomplished by: 
(1) altering program statements necessary to reconcile the difference in 
word lengths between the two types of computers, (2) making various 
syntax changes necessitated by technical differences in the implementa- 
tion of FORTRAN on the different machines, and (3) changing the modular 
form and adding program statements necessary to provide for the execution 
of the algorithm as intended by the author. Details of the conversion are 
apparent from comparison of the program listing available at the W. R. 
Church Computer Center, USNPGS, Monterey, California with that of Ref. 20. 
The listing found in the appendix of this thesis contains only the NAS], 
NAS4, and SENS subroutines. The first two are included since extensive 
modifications were incorporated in the course of implementing the 
sensitivity feature, and the latter because it is entirely new. 

As previously described, the NASAP program has available, in symbolic 
form, the uncancelled loop products created from a topological flowgraph 
analysis of a given circuit. It was from this base that the writer 
devised a program to compute circuit transfer function sensitivities in 
an algorithm developed as follows. 

Since, in subroutine NAS4 of the original program, there implicitly 
exists a requested transfer function, 


T(s) = K N(s) 
D(s 





in which, 
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Gmere nN, M, p, q, n',m', p', q', ... are integer exponents of the Laplace 


variable s, not sorted into any ascending or descending order and with 


th 


repetitions possible; and where as» b5 represent the i loop products 


(algebraic product/quotients of circuit parameters in symbolic form); 
and where LPD is the number of uncancelled loop products in the denomi- 
nator polynomial of T(s) and LPN the number in the numerator. Each loop 
product is actually determined from one 32-bit computer word, wherein 
the presence (1) or absence (0) of each element is indicated by the 
corresponding bit value. This, of course limits the number of circuit 
parameters which can be accommodated by the program to NBITS-1 = 31, 
although in the current version this upper limit is set to 23 due to other 
considerations. | 

It was desired, then, to compute the sensitivity of T(s) in symbolic 


form, to a specified circuit parameter P Since sensitivity is essen- 


kK! 


tially a rate-of-change, it follows that the transfer function sensitivity 





may be “ as: | D(s) .oNCs) a N(s) 22{s) 
(gs, = aM(s) = a(N(s)/D(s)) - ‘ 
3P 3 P ; 2 (1) 
k k k D (s) 


To implement this fundamental approach utilizing the data available 


in NAS4, the following simple algorithm was developed, (See Flow-Chart, 
Eee. 1): 

_th th 
(I) Observe the j factor of the i loop product, b,, , in the numerator 
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(II) Repeat step (I) for the Genet factor of the ith loop product. 
Continue until all factors are exhausted (j>NFN(i), where NFN(i) 


represents the number of factors in the sie loop product of the 


numerator polynomial). 


(III) Does es contain PL? (is flag set tol?) 
Kk 


YES: Continue 
NO: Set an | =e . 
oPy ; 
i 


IV) (A) i=itl 
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(B) Repeat (I) (II), and (III) for the (i41)t® loop product. 


Continue until all loops are exhausted (i>LPN). 
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(V) Repeat the entire procedure beginning with step (I) for the 
deneminator of Tis). 
At this point the sensitivity is implicitly available in the desired form 


(eqn. |), where: 
LPD NFD(i) EXPDS (i) 
D Ss 


D = » II 4 
i=l j=l zs 
LPN NFN(i) EXPNS (i) 
N _ suds JT wpa 
i=l j=l J 
LPD NFD(i) ji 
Ey - 2 0 (= | _EXPDS(i) 
1= = 
e J 3Py ne 


NEN (3 E t 
ON ue - Gy (2 \ 5 EXPNS (4) 
ij 


oP i= — 
k i=l] jel oP, 


h 
NFN(i) = Number of factors in the i* loop product of numerator 


polynomial 
NFD(i) = Number of factors in the gee loop product of denominator 
EXPNS(i) = Exponent of s in the th loop product of numerator 
EXPDS(i) = Exponent of s in the zth loop product of denominator 


One feature of the programming which became necessary in order to 
reduce storage requirements for the symbolic sensitivity component variables 
SNUM, SDEN, PARTLN, and PARTLD was a packing scheme which in effect 
reduced vector sizes by a factor of one-fourth. In the NASAP program, as 
in most FORTRAN programs, each alphanumeric character is represented by 
a four-byte word, wasting a great deal of core storage since one byte is 
sufficient for one character. This exemplifies one of the most unattractive 
aspects of the FORTRAN language, especially when dealing with character 
Strings or literals as in the case of the NASAP symbolic version. 

A remedy for the problem which was utilized in this program is a 


very short but efficient packing subroutine which has four arguments: 
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(1) a four-byte variable to be "packed" 
(2) an integer constant specifying the byte number of the variable 
in (1) to be assigned the value of the packing variable 
(3) packing variable 
(4) an integer constant specifying the byte number of the variable 
in (3) which is to be assigned to the variable in (1) 
Upon calling the pack subroutine, dummy variables, locally declared 
LOGICAL*1 (one byte), receive data passed from the specified byte of the 
packing variable, assign this to the specified byte of the variable to be 
packed, and then return this new "packed" version to the calling program. 
In this manner the sensitivity component variables were packed with para- 
meter name (e.g. LO1, RO4, etc.) in bytes one through three, and byte 
number four was packed with a code number to indicate the form in which 
the parameter should appear. This code appears below: 


O -- Numerator or denominator 


= 7Z, partial = z 
3 -- Numerator or denominator = z, partial = 1 
4 -- Numerator or denominator = z, partial = 0 
6 ~~ Numerator or denominator = 1/z, partial = 1/z 


7 -- Numerator or denominator 


1/z, partial = -1/22 


l/z, pastial = 0 


9 -- Numerator or denominator 

With this scheme the program was rendered capable of processing net- 
works of up to 200 independent loops to provide for-sensitivity with 
respect to a single parameter, or by successive solutions to provide for 
any or all sensitivities desired. With this loop-set limit, the entire 
program, including NASAP requires 180K of core. A guide fox core storage 


allocation may be found in the appendix for circuits of greater complexity. 


14 





Subsequent to the packing operation, data is transferred to a new 
Subroutine SENS for serait on into the desired form and for numerical 
evaluation. In this operation it was found impractical, of course, to 
permit each product element to occupy a separate storage location in the 
computer, Since this would essentially require dimensions on the order 
of the square of the maximum number of loop products. Therefore, numerical 
evaluation of the constituent components of the sensitivity are recursively 
computed and simultaneously combined at this point in the program to form 
the numerator and denominator arrays, PART1 and PART2. Eliminating all 
zero vectors and combining like powers of s reduced storage requirements 
for this operation essentially to twice the order of the denominator | 


polynomial. 
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IIT. USER'S RULES 


Users of the sensitivity subroutine must be familiar with the simple 
NASAP coding rules contained in Refs. 1] and 2 as modified below: 

(1) Only one "OUTPUT" and/or "ROOTS" transfer function request should 
be made in a given problem since the sensitivity subroutine is 
designed to store data relevant to only one numerator polynomial 
at a time. Of course, the user may wish to follow one set of 
problem cards with another. 

(2) The "SYMBOLIC SOLUTION" card must be inserted following "EXECUTE". 

(3) No "MODIFY" requests are permitted. 

(4) The card following USYMBOLLC SOLUTION must cottain the words 
(only the underlined portions are essential for machine 
recognition) ; 

SENSITIVITY PARAMETER=LO2 (or RO1, CO9,. etc.) 
Any number of blanks or identifying characters may be included 
so long as the underlined portions appear somewhere in columns 
1-80 in the order shown. Notice that the sacemate: (with respect 
to which the sensitivity is desired) is indicated by the consecu- 
tive characters, (=), (R, L, or C), (two digit element number). 
(5) If it is desired to have the sensitivity magnitude and phase 
evaluated, the next card in sequence must contain a maximum 
frequency in columns 1-11 and a minimum frequency (not zero) in 
columns 21-31, preferably both in standard FORTRAN formal E11.4, 


although normal floating-point form is satisfactory if all 


io 





significant digits (plus sign) remain within the specified 
11 columns. Specifying these frequencies limits will result in a 
print-out of the sensitivity magnitude and phase at 200 evenly 
distributed frequencies in the interval. 
(Caution: If the user omits this option and desires to follows 
the current problem deck with others, a blank card must be 
inserted in place of the frequency interval card). 
An interested user may provide for the plotting of the generated 

magnitude and phase data by changing the variables RMAGN, ANGLE, and 

FREQ into array variables in the SENS subroutine, then passing them to 

an existing library plotting subroutine, This will, of course, consider- 


ably increase the core storage required for the program. 
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IV. INTERACTIVE APPLICATIONS 


The NASAP program, in conjunction with the sensitivity subroutine, 
may be easily adapted for use under a time-sharing system such as CP/CMS 
with only a few simple statement changes which depend upon the user's 
choice of file structure. One choice successfully used at the Naval 
Postgraduate School is to remove the four initial statements: (1) CALL 
NASAP, (2) STOP, (3) END, and (4) SUBROUTINE NASAP. It is possible, then, 


to "offline read" the program into CP in the following modular form: 


FILE PROGRAM CONTENT 

RSS Program RSS (Driving program) 
NAS1 NAS1, FUNCTION ITRAN 

NAS2 NAS2 

NAS3 NAS 3 

NAS4 NAS4, SUBROUTINE PACK 

NAS5 NASS5 

SENS SENS, BLOCK DATA 


Under CP/CMS, the entire program can then be implemented by executing 
Miss, with data being input at the terminal, or if the three "READ" 
statements in NASI are changed appropriately, data may be read from a 


separate file so defined. 


al 





V. DISCUSSION AND RECOMMENDATIONS FOR FURTHER IMPROVEMENT 


As the reader may have gathered, the feasibility of computing sensi- 
tivities of a transfer function to all circuit parameters and storing 
them in symbolic form appears doubtful under present storage limitations. 
What follows, however, is a recommendation for further reducing storage 
requirements somewhat while maintaining the symbolic concept. If only 
a single numerical sensitivity were desired, it would of course be 
considerably more efficient to recursively evaluate the sensitivity 
from successive loop-sets with their appropriate values first as given, 
then in partial derivative form. This procedure would eliminate most 
of the storage problems. The singular advantage of the symbolic approach, 
however, is that evaluation of the sensitivity function at various fre- 
quencies, of for various parameters, is a very simple matter of substi- 
tution, rather than the use of time-consuming iterative or numerical 
techniques which may pose stability or accuracy difficulties. 

It is apparent that substantial changes are indicated in the proposed 
algorithm (or perhaps only the techniques for implementing the algorithm 
which will permit the sensitivity subroutine to cope with circuits of 
greater complexity by means of a reduction in the large storage require- 
ments. However, it must be remembered that the resulting exponential 
increase in the number of independent loops with respect to a concomi- 
tant increase in the number of nodes [20] will in itself limit the 
complexity of topological/flowgraph methods, particularly in maximal- 
planar networks. Indeed, the difficulties of dimensionality plague all 


circuit analysis programs. 
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One suggested method of reducing storage requirements would be 
to extend the previously mentioned concept of utilizing individual 
bits of the computer word to represent constituent elements of the 
network. Since in this case it is necessary to know which of six 
mathematical formulations is appropriate for each element of the 
numerator and denominator polynomials, instead of only two formula- 
tions as in the case of transfer function loops, it would appear 
that a scheme requiring only three variables (23 = 8 possibilities), for 
a numerator and denominator each dimensioned to the maximum number of 
loops (10,000 in the present version) would suffice to allow sensitivi- 
ties to be computed and stored in symbolic form. This would repre- 
sent a substantial reduction from the 460,000 that are presently required 
for a 10,000-loop circuit to accomplish the same task. This suggestion, 
in effect, offers an execution time-vs-storage tradeoff, since a coding/ 
decoding algorithm must be developed which will intrinsically involve 
a tremendous number of manipulations. In addition it would require an 
almost total rewrite of the FORTRAN programming, particularly that 
portion meshed with the NASAP NAS4 subroutine, and it may very well 


lead to difficulties not apparent without further investigation. 
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APPENDIX A SAMPLE PROBLEMS 


(I) Sample problem No. 1 - Emitter Follower 


(a) Emitter follower linear model: 





(b) NASAP coding: 


NASAP PROBLEM NO. 1 
xEMITTER FOLLOWER 


v1121 

a 1 ow 1 

pa 2 3 M00 
R23 4 1M 
R33 11.5K 
Ru & 1 100 
Cl 3 4 6PF 
C2 3 1:200PF 
12 & 1 0.04VC2 
OUTPUT 
VRU/VV1 
EXECUTE 


SYMBOLIC SOLUTION 
SENSITIVITY PARAMETER=R0O4 
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(c) Sensitivity output for emitter follower: 


NUMERATOR POLYNOMIAL : SENSITIVITY OF VRO4/VVO1 TO ROK 
O.6944E 40 TIMES S** 0 
0.1042E 31 TIMES S** -1 


0.6667E 17 TIMES S** -2 


DENOMINATOR POLYNOMIAL 3; SENSITIVITY 


Oey se 19 TIMES S**. 2 
O54 75E Sat IMES Seo 2] 


QonhvesE HS -IRIES S** 0 


O.1111E 30 TIMES S** -] 
0.2845E 16 TIMES S** -2 
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(II) Sample Problem No. 2 - Butterworth Filter 


(a) Original linear circuit: 





(b) NASAP coding: 


NASAP PROBLEM 
*THIRD ORDER BUTTERWORTH 
R1 4 2 ] 


i 
1 
2 
| 
‘| 


SYMBOLIC SOLUTION 
SENSITIVITY PARAMETER=R02 
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(c) Sensitivity output for Sample Problem No. 2 - Butterworth 


Filter: 


NUMERATOR POLYNOMIAL : SENSITIVITY OF VRO2/VVO1 TO ROQ2 
Q.5000E 00 TIMES S** 0 
0.5000E O00 TIMES S** -] 


0.2500E 00 TIMES S** .-2 


DENOMINATOR POLYNOMIAL : SENSITIVITY 
O,1000E O1 TIMES S** 4& 

OVHUOUE OL TBs ”S** 35 

0.8000E O1 TIMES S** 2 

LTUOOE 02 THES S** 1 

O0.8000E O1 TIMES S** 0 


O,.4000E O1 TIMES S** -] 
De tOUgEsO] TIMES S** =2 


ae | 





(III) Sample Problem No. 


(a) Original circuit 


3 —- Transformer circuit 


with ideal transformer: 





(b) Linear model: 





(c) NASAP coding: 


NASAP PROBLEM 


*|NPUT 
fie 2 
R12 4 
GL > 1 
fi 2) > 
[2 oa i 
Wee He 
OUTPUT 


VIi/ttl 
EXECUTE 


IMPEDANCE 
1 

> 

0.5 

0.25 

=e eb 

=2 VEL 
REGUL S fi 


SYMBOLIC SOLUTION 
SENSITIVITY PARAMETER=LO1 
woe 05 ~1592E 00 
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(d) Sensitivity output for Sample problem No. 3 - Transformer 


CireuLt: 


NUMERATOR POLYNOMIAL : SENSITIVITY OF VIOI/!I!01 TO LO] 
OneeOOOUE OF TIMES S** |] 
~Oeoeoor 00 TRIES S*= ©@ 


0.2778E-01 TIMES S** -1 


DENOMINATOR POLYNOMIAL : SENSITIVITY 
‘OR lyee OF TIMES Swe 2 

0.2667E O1 TIMES S** 1] 

Opiate OL TRMES Seer 0 

OFSsooke O00 TRMES Se -1 


0.2778E-01 TIMES S** -2 
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(e) 


Sensitivity magnitude and phase output for Sample Problem 


No. 3 - Transformer circuit: 


EVALUATION OF SENSITIVITY 


FREQUENCY 


0,1592E 
0.9584hE 
0. 1L758E 
0.255 7E 
0.3356E 
OQ) , 4 SIE 
O.4954E 
ORG Se 
Oe 6555 
02-755 25 
0. S151e 
'0.8950E 


00 
00 
01 
Ol 
Ol 
jal 
01 
O01 
O01 
O01 
O01 
Ol 


OU.O2Z5U0E 
O65 50'E 
0.6409E 
0.6489E 
O65 69E 
0.6649E 
w.67Z9E 
0. boug'eE 
10.6889E 
O36 9 JE 


0.1480E 
0.1488E 
0.1496E 
O.1504E 
0.1512E 
0.1520E 
0.15286 
0.1536E 
O.1544E 
0.1552E 
0.1560E 
0.1568E 
0.1576E 
0.1584E 
0.1592E 


02 
02 
02 
02 
02 
02 
02 
02 
02 
02 


MAGNITUDE 


OW. 4552E 00 
Om Zb8E> 01 
Ue oUe7ZeE-Ol 
Deets 8SE= 01 
0, 2666E-01 
Weal te= 0) 
Oe ovibE= O01 
Ue oi = 01 
UF labo e= 01 
DO. PABsE-01 
0, wu ok 01 
0.1000E-01 


0.1432E-02 
0.1414E-02 
0.1397E-02 
0:1380E-02 
0.1363E-02 
0.1346E-02 
0.1330E-02 
0.1315E-02 
0.1300E-02 
0.1285E-02 


0.6049E-03 
0.6016E-03 
0.5984E-03 
0.5952E-03 
0.5921E-03 
0.5890E-03 
0.5859E-03 
0.5828E-03 
0.5798E-03 
0.5768E-03 
0.5739E-03 
0.5709E-03 
0.5680E-03 
0.5652E-03 
0.5623E-03 
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eed 


PHASE CIN 


ORE) bE 
=Oiely2.b / E 
po gelled ts OS 
~O2GS7E 
-O.,1484E 
=O77 POLE 
=0 son 2 E 
=) Ser OE 
=O Po Z0c 
, OR OL 
it. DE 
=). OE 


a 


=0 ebb E 


=D D OWE 
=O, 15666 
-0..47 > 668 


-0.1566E 
-0.1566E 
-0.1566E 
-0.1567E 


-0.1567E 
-0.1567E 


=). LoOge 
oO) edt) 6) Sie 
ee) Salles eek 
Osu o ot 
0% bo 65 
ge lg desi 6 818 
ely de slaye)(e 
=O, 15636 
=O), 5 691 
+00 © Se 
Oreo 6 he 
=O aol GE 
=o 
SO or Je 
pa Uh ee os 


RADIANS) 


00 
Ol 
O01 





(IV) Sample Problem No. 4 - I.F. Transformer 


Cy &£.Fo Transformer circuit: 





(b) Linear model: 





(c) NASAP coding: 


NASAP PROBLEM 
*[.F. TRANSFORMER 


falas 2 oe 

Ri 2°). 
Cie 7a) 2S Pe 
Rize 2 B10 

Li 3 & .4&9MH 
L2 4 1 .O11H 
L3 4& 5 .4&9MH 
R3 5 6 10 

C2 6 1 2h0PF 
Gor 2° 6°2 PF 
R4 6 1 15M 
OUTPUT 

VR4E/ TI 
EXECUTE 


SYMBOLIC SOLUTION: 
SENSITIVITY PARAMETER=C03 
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(d) 


sensitivity output for Sample Problem No. 4 - I.F. 


Transformer: 


- NUMERATOR POLYNOMIAL 


0.4720E 
0.4786E 
0.1293E 
0.4049E 
0.4725E 
0.2309E 


Oe 4 VUGE- 


0.4760E 


WLS 2 7 E 


0.1304E 
0.4095E 
0.4805E 
0.2375E 
0.4266E 
0.4074E 


ORL OOCE 


18 
16 
14 
10 
06 
02 
05 


DENOMINATOR 


42 
40 
38 
34 
30 
Z6 
7 J. 
15 


09 


TIMES 
Ties 
TIMES 
TIMES 
TIMES 
TIMES 
TIMES 


POLYNOMIAL 


TIMES 
TIMES 
TIMES 
eS 
TIMES 
TIMES 
TIMES 
TIMES 


TIMES 


S«** 


S** 


S** 


S** 


S** 


S*e 


S** 


S** 


S** 


S** 


S*x* 


S** 


S** 


S** 


S*x* 


S** 


SENSITIVITY OF VROS/TIOL 


> 


6 


5 


4 
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TO 
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(V) Sample Problem No. 5 - Wien-Bridge circuit 


()Oraeinal circuit: 


. 


(b) Linear model: . oan 
con 






(c) GNASAP scoding: 


NASAP PROBLE}M 

‘VW EN-BRIDGE CIRCUIT 
ile 2 
ale 2 
Caz 
Ree 
C25 > 
N Goel 
OUNGeU T 

VV2/111 

ESeCUTE 

SYMBOLIC SOLUTION 


SENSITIVITY PARAMETER=R02 


— FN pk 
i bt ke pe OF 
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(d) Sensitivity output for Sample Problem No. 5 - Wien-Bridge 


cirecult: 


ROMERATOR  POLYNOMTAL : SENSITIVITY OF VVO02/71101 TO RO2 


IDENTICALLY ZERO 


DENOMINATOR POLYNOMEAL : SENSITIVITY 


OFEOO0E O1 TIMERS Sa xggek 


O.4000F O01 TIMES Sx 
U,G000F O1 TIMES S** 


Oyremoe O1 TIMES S** 1 


OeeOUoE Ol Tits s*a 0 
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(V1) SanpiiesProbliem No. 6 - High-Q Filter 


Coop riiter circuit: 





(b) NASAP Coding: 


NASAP PROBLEM 
mine sie WE 


Mellor eel 

R1 2 4 20K 

RZ Oe” Deal OK, 

R3 6 & 20K 

R4 6 7 100 

Ro 7 5 WOO 

R6 5 1 2K 

Ry 7 1 Ti 

C12 Saaeto 2 
C2 5 tO UE 
CanG 36.0 2UF 
er le? low tans 
ieee 5150 IRS 
Clrr rw T 

VR6/VV1 
EXECUTE 


SYMBOLIC SOLUTION 
SENSITIVITY PARAMETER=C02 
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(c) Sensitivity output for Sample Problem No. 6 — High-Q filter: 


NUMERATOR POLYNOMIAL : SENSITIVITY OF VROG/VVO01 TO CO2 
~0.2684E 09 TIMES S** 6 | | 
Ueoteor el? TIMES S** 5 
O,S125— 09 TIMES Se* 4 
O.1031E 06 TIMES S** 3 
fOeeleote O2 TUMES S#&°" 2 


.0.5000E-03 TIMES S** 1 


DENOMINATOR POLYNOMIAL : SENSITIVITY 
0.1650£ 24 TIMES S** 6 
0.2869E 21 TIMES S** 5 
0.1887E 18 TIMES S** 4 
0.5643E 14 TIMES S** 3 
0.6908E 10 TIMES S** 2 
0.1575E 06 TIMES S** 1 
0.1000F 01 TIMES S** 0 
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APPENDIX B GUIDE FOR THE ALLOCATION OF CORE STORAGE 


NUMBER OF LOOP-SETS POSSIBLE 


(VARIABLE DIMENSION) PROGRAM CORE REQUIRED 
* 200 180K 
300 205K 
400 230K 
500 255K 
600 280K 
700 305K 
** 800 330K 
900 355K 
1000 380K 
*kk 5000 505K 
10000 630K 
* Corresponds approximately to a maximal planar network 


of 5 nodes 

KK Corresponds approximately to a maximal planar network 
of 6 nodes 

ARK Corresponds approximately to a maximal planar network 


of 7 nodes 


ay 





PROGRAM LISTING 


FRMIN 
),ITS 


VY) we SOY 
as ofp <i Y 
(Ace SS” 
Ivy J<Ir: 
2 Oo. “. 
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rom NY ow 
DWNZZrROZ= 
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prio L023 


I 
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1361 


Way TSO) + ING (3b e422) > Nanna 
ea) Te 


PNEST 
Bee and 2 


=1£381300 514300 
Q.1SYMB(MM) ) 


298 } 
YMB(46))1353,1355,1353 


7 a1 E/E 


ND 


793032352168 168)sNASA 
3-80 
15 


01,351);1 


ODOR R ae rw 
wes ONY) ION DAFNRHONDANNDM eODtHO 
Vo ae See ae 1 oe ee Sg We Wo US Del © Ds De] | a | 
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eve eee le ae ee ee ee ee 
fod mt OH ee ee 


Ns OR ER SE 
MOO -INWnNND 
Tar ina ea 


erage a 


300 
250 
P26 


I. 
1,10 


1+2}.£Q.1SYMB(MM)) GO TO 1363 


M 
M 
MM-1 


ta 


Or 
0 0 
mr) 
onKamn! 


FRMAX,FRMIN 
38 


60) 
68;68169)+NASA 


313051 68:68)sNASA 


O, END 
159% 
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